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ABSTRACT 

Porcine deltacoronavirus (PDCoV) is an emerging swine enteropathogenic 
coronavirus. The first outbreak of PDCoV was announced from the United States in 
2014, followed by reports in Asia. The nonstructural protein nsp5 is a 3C-like 
protease of coronavirus and our previous study showed that PDCoV nsp5 inhibits type 
I interferon (IFN) production. In this study, we found that PDCoV nsp5 significantly 
inhibited IFN-stimulated response element (ISRE) promoter activity and transcription 
of IFN-stimulated genes (ISGs), suggesting that PDCoV nsp5 also suppresses IFN 
signaling. Detailed analysis showed that nsp5 cleaved signal transducer and activator 
of transcription 2 (STAT2), but not Janus kinase 1 (JAK1), tyrosine kinase 2 (TYK2), 
STATI and interferon regulatory factor 9 (IRF9), key molecules of the JAK-STAT 
pathway. STAT2 cleavage was dependent on the protease activity of nsp5. 
Interestingly, nsp5 cleaved STAT2 at two sites, glutamine (Q) 685 and Q758, and 
similar cleavage was observed in PDCoV-infected cells. As expected, cleaved STAT2 
impaired the ability to induce ISGs, demonstrating that STAT2 cleavage is an 
important mechanism utilized by PDCoV nsp5 to antagonize IFN signaling. We also 
discussed the substrate selection and binding mode of PDCoV nsp5 by homologous 
modeling of PDCoV nsp5 with the two cleaved peptide substrates. Taken together, our 
study demonstrates that PDCoV nsp5 antagonizes type I IFN signaling by cleaving 
STAT2 and provides structural insights to comprehend the cleavage mechanism of 
PDCoV nsp5, revealing a potential new function for PDCoV nsp5 in type I IFN 


signaling. 
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IMPORTANCE 

The 3C-like protease encoded by nsp5 is a major protease of coronaviruses; thus it is 
an attractive target for development of anti-coronavirus drugs. Previous studies have 
revealed that the 3C-like protease of coronaviruses, including PDCoV and porcine 
epidemic diarrhea virus (PEDV), antagonizes type I IFN production by targeting 
NF-«B essential modulator (NEMO). Here, for the first time, we demonstrate that 
overexpression of PDCoV nsp5 also antagonizes IFN signaling by cleaving STAT2, 
an essential component of transcription factor complex ISGF3, and that PDCoV 


infection reduces the levels of STAT2, which may affect the innate immune response. 
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INTRODUCTION 

Porcine deltacoronavirus (PDCoV) is an emerging swine enteropathogenic 
coronavirus (CoV), which causes acute diarrhea, vomiting, and mortality in piglets 
(1-4). PDCoV was first found from pig feces in Hong Kong during molecular 
surveillance of CoVs in avian and mammalian species in 2012 (5). The first outbreak 
of PDCoV was announced in the United States in 2014 (6). To date, it has been 
detected in at least 20 states of the United States (3, 7, 8), as well as Canada, South 
Korea (9), China (10-12), Thailand (13), Lao People's Democratic Republic, and 
Vietnam (14), posing a significant threat to the swine industry and gaining 
considerable attention (13, 15, 16). 

PDCoV is an enveloped, single-stranded positive-sense RNA virus and belongs 
to the genus Deltacoronavirus in the family Coronaviridae (5, 17). The PDCoV 
genome, approximately 25 kb in length, encodes four structural proteins, spike (S), 
envelope (E), membrane (M), and nucleocapsid (N), and three accessory proteins, 
NS6, NS7, and NS7a, interspersed within the structural genes at the 3’ end of the 
genome (5, 18, 19). Open reading frame (ORF)1a and ORF 1b are located downstream 
from the 5’ untranslated region and encode two polyprotein precursors that are 
predicted to be processed into 15 mature nonstructural proteins (nsps) by the 
papain-like protease, nsp3 and the 3C-like protease, nsp5 (20, 21). Because of the 
essential function of nsp5 in viral replication, it has become an attractive target for the 
development of anti-CoV drugs (22, 23). Understanding in detail the function and 


structural characteristics of PDCoV nsp5 will aid in development of targeted therapies 
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against CoV infection. 

Interferon (IFN) and the IFN-induced cellular antiviral response are primary 
defense mechanisms against viral infection (24). Type I IFN binds to its surface 
receptor at two subunits, IFNARI and IFNAR2, and induces heterodimerization of 
the receptor subunits. Subsequently, Janus kinase 1 (JAK1) and tyrosine kinase 2 
(TYK2) are activated and phosphorylate the cytoplasmic domains of IFNAR1 and 
IFNAR2? to establish docking sites for signal transducer and activator of transcription 
1 (STATI) and STAT2 (25, 26). STAT1 and STAT2 are then phosphorylated by JAK1 
and TYK2, leading to heterodimer formation of phosphorylated STAT! and 
phosphorylated STAT2 (26, 27). This complex interacts with IFN regulatory factor 9 
(IRF9) to form IFN-stimulated gene factor 3 (ISGF3), which is transported into the 
nucleus and recognizes IFN-stimulated response elements (ISREs), enhancing the 
transcription of hundreds of IFN-stimulated genes (ISGs) (28, 29). As effective 
antiviral factors, ISGs directly lead to resistance to viral infection by targeting viral 
processes such as transcription and translation. 

Given the importance of JAK-STAT pathway in the IFN signaling, it is not 
surprising that viruses have evolved mechanisms that target it to inhibit ISG 
production. Many studies have demonstrated that viral proteins, such as the NS5 
protein of dengue virus (DENV), Zika virus (ZIKV), and the hepatitis C virus (HCV) 
core protein, target the JAK-STAT pathway by diverse mechanisms to negatively 
modulate IFN signaling (30-32). We previously demonstrated that PDCoV nsp5 


antagonizes IFN-B production by cleaving NF-«B essential modulator (NEMO) (33). 
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However, whether nsp5 regulates IFN signaling remains unclear. In this study, we 
found that PDCoV nsp5 suppressed type I IFN signaling by cleaving STAT2, a key 
molecule in the JAK-STAT pathway, revealing a potential new function for PDCoV 


nsp5 in type I IFN signaling. 


RESULTS 

Identification of PDCoV nsp5 as an antagonist of type I IFN signaling. Type I 
IFN signaling induces a potent antiviral state by enhancing the expression of hundreds 
of ISGs, which are critical for controlling viral infections (34). To assess the 
regulatory role of PDCoV nsp5 in type I IFN signaling, mRNA levels of ISG15, 
ISG54, ISG56, and OAS1 were evaluated under PDCoV nsp5 overexpression in 
human embryonic kidney cells (HEK-293T). As shown in Fig. 1A, ISG transcription 
induced by IFN-a was significantly inhibited by PDCoV nsp5 compared with the 
control group. Because of the existence of ISREs within ISG promoter regions, 
various concentrations of PDCoV nsp5 expression plasmid and ISRE-luciferase 
reporter plasmid were co-transfected into HEK-293T cells, porcine kidney cells 
(PK-15) or swine testicular cells (ST). The results showed that nsp5 strongly inhibited 
IFN-a-induced ISRE promoter activity in a dose-dependent manner in HEK-293T 
cells (Fig. 1B), PK-15 cells (Fig. 1C) and ST cells (Fig. 1D). These results confirm 


the antagonistic property of PDCoV nsp5 in type I IFN signaling. 


PDCoV nsp5-mediated inhibition of type I IFN signaling requires its protease 
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activity. Because cysteine (Cys) and histidine (His) residues of CoV nsp5 form a 
catalytic dyad, any mutation in the catalytic site can disrupt its protease activity (22, 
35, 36). Secondary structure alignment showed that residues His41 and Cys144 
(numbering based on PDCoV nsp5) are highly conserved among the four genera of 
alpha-, beta-, gamma-, and delta-CoVs (Fig. 2A). We further created a homology 
model of PDCoV nsp5 based on the crystal structures of porcine epidemic diarrhea 
virus (PEDV) nsp5 (PDB identifier 4ZUH) using MODELLER 9.14 software and 
compared it with the three-dimensional (3D) structures of PEDV nsp5 (alpha-CoV), 
severe acute respiratory syndrome coronavirus (SARS-CoV) nsp5 (beta-CoV), and 
infectious bronchitis virus (IBV) nsp5 (gamma-CoV). As shown in Fig. 2B, the 
locations of the catalytic residues are also conserved in the 3D structures of PEDV 
nsp5, SARS-CoV nsp5, IBV nsp5, and PDCoV nsp5. Thus, two mutations, nsp5 
H41A and nsp5 C144A, in the catalytic dyad were constructed. Compared with 
overexpression of wild-type PDCoV nsp5, repression of IFN-a-induced ISRE 
promoter activity was strongly relieved with overexpression of nsp5 H41A or nsp5 
C144A (Fig. 2C). This finding suggests that the protease activity of PDCoV nsp5 is 


involved in antagonizing type I IFN signaling. 


PDCoV nsp5 antagonizes type I IFN signaling by targeting STAT2. To determine 
the mechanism by which PDCoV nsp5 inhibits type I IFN signaling, the endogenous 


protein level and phosphorylation of JAK1, TYK2, STATI, and STAT2 were 


examined in IFN-a-treated LLC-PK1 cells (porcine kidney cell line). As shown in Fig. 
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3A, the expression and phosphorylation of JAK1 and TYK2 were not affected by 
PDCoV nsp5. However, STAT2 phosphorylation was impaired by PDCoV nsp5S after 
IFN-a treatment and the expression of total STAT2 was reduced (Fig. 3B). Thus, we 
did not exclude the possibility that the weaker band of P-STAT2 in the presence of 
nsp5 is probably due to the reduced levels of total STAT2. In addition, a slight 
reduction in endogenous STAT] expression was also observed (Fig. 3B). 

As for the importance of PDCoV nsp5 protease activity in antagonizing type I 
IFN signaling, we considered whether PDCoV nsp5 cleaves specific molecules. 
Therefore, nsp5 was transfected into HEK-293T cells, along with JAK1, TYK2, 
STATI, STAT2, or IRF9. No cleaved products were detected with JAK1, TYK2, 
STATI, or IRF9 cotransfection. However, two faster migrating protein bands were 
observed by western blot analysis in cells cotransfected with STAT2 and PDCoV nspS, 
indicating that PDCoV nsp5 mediates cleavage of STAT2 (Fig. 3C). Although nsp5 
overexpression also slightly reduced the protein abundance of endogenous STAT1 
(Fig. 3B), no cleavage product of Flag-STAT1 in nsp5-transfected cells was detected 
by an anti-Flag antibody (Fig. 3C). The nsp5-mediated reduction in endogenous 


STATI expression was not investigated further in this study. 


PDCoV nsp5 cleaves STAT2 through its protease activity. To further confirm nsp5 
cleavage of STAT2, HEK-293T cells were cotransfected with STAT2 and various 
concentrations of PDCoV nsp5. PDCoV nsp5 induced STAT2 cleavage in a 


dose-dependent manner (Fig. 4A). Additionally, STAT2 cleavage was not observed in 
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cells cotransfected with nsp5 H41A or nsp5 C144A (Fig. 4B). Furthermore, PDCoV 
nsp5-induced STAT2 cleavage was not affected by treatment with caspase inhibitor 
Z-VAD-FMK or proteasome inhibitor MG132, confirming this cleavage occurred 
independent of cellular caspases or the proteasome (Fig. 4C). Taken together, PDCoV 
nsp5 cleavage of STAT2 appears to occur through its protease activity. 

To further evaluate the expression of endogenous STAT2 in PDCoV infection, 
the mRNA and protein levels of endogenous STAT2 were determined in 
PDCoV-infected LLC-PK1 cells. Although STAT2 mRNA in PDCoV-infected cells 
slightly increased compared with the mock-infected group (Fig. 4D), endogenous 
STAT2 protein was clearly decreased with PDCoV infection (Fig. 4E). Remarkably, 
the cleavage of endogenous STAT2 was also detected in PDCoV-infected cells (Fig. 


4E), revealing the biological effect of this cleavage in virus infection. 


PDCoV nsp5 cleaves STAT2 at residues Q685 and Q758. To further identify the 
cleavage sites of STAT2 by PDCoV nsp5, substrate specificity of PDCoV nsp5 in 
viral polyprotein itself was analyzed. As shown in Fig. 5A, a preference for substrate 
cleavage by PDCoV nsp5 was Q residue at the P1 position, which is common among 
other CoV subfamilies (1, 37). Because the cleavage of N-terminally Flag-tagged 
STAT2 by PDCoV nsp5 yielded a slightly longer than 70 kDa in Fig. 3C detected by 
anti-Flag antibody, a cluster of potential nsp5 cleavage sites were suggested to be in 
the carboxyl-terminal region of STAT2. Considering the molecule weight of cleavage 


products and the conserved Q residue at the Pl position in substrate recognition by 
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coronavirus nsp5, the deletion mutants, STAT2 (1-637 aa), STAT2 (1-719 aa), and 
STAT2 (1-758 aa) were generated. The two cleavage products of STAT2 were longer 
in size than STAT2 (1-637 aa), indicating that the recognized Q residue of the P1 
position should be located at the C-terminal region (residues 637-864 aa) of STAT2 
(Fig. 5B). Eight Q residues (637—864 aa region) were subsequently changed to alanine 
(Fig. 5C) and transfected into HEK-293T cells along with PDCoV nsp5. The two 
cleavage products separately disappeared under the expression of Q685A or Q758A 


(Fig. 5D), while other mutants were cleaved into two new products as expected by 


PDCoV nsp5, demonstrating the two potential cleaved positions to be Q685 and Q758. 


To further confirm these cleaved residues, a double mutant (STAT2-Q685A-Q758A) 
was constructed and cotransfected into HEK-293T cells with PDCoV nsp5. The 
results showed that STAT2-Q685A-Q758A was resistant to cleavage, confirming 
Q685 and Q758 to be the two P1 position residues recognized by PDCoV nsp5 for 


cleavage of STAT2 (Fig. 5E). 


PDCoV nsp5-mediated cleavage impairs STAT2’s ability to induce ISGs. To 
assess the effect of STAT2 cleavage by PDCoV nsp5 in type I IFN signaling, several 
deletion mutations encoding STAT2 cleavage products, including STAT2 (1-685 aa), 
STAT2 (1-758 aa), STAT2 (686-864 aa), STAT2 (759-864 aa), and STAT2 (686-758 
aa), were generated based on the identified cleaved residues. Full-length STAT2, 
individual deletion mutants, or mixed deletion mutants were cotransfected with 
STAT1 and IRF9, along with ISRE-luciferase reporter plasmid and pRL-TK. 
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Compared with full-length STAT2, the luciferase activity of ISRE-Luc significantly 
decreased under the expression of any cleaved fragments (Fig. 6A). The mRNA levels 
of ISGs, such as ISG15 (Fig. 6B), ISG54 (Fig. 6C) and OAS1 (Fig. 6D), were also not 
sufficiently induced by the cleaved STAT2 fragments. These results reveal that STAT2 
cleavage induced by PDCoV nsp5 impairs the ability of STAT2 to activate type I IFN 


signaling. 


DISCUSSION 

As a powerful antiviral response, type I IFN signaling can control viral infections 
through activation of the transcription factor complex ISGF3, resulting in increased 
transcription of hundreds of ISGs that help to generate an antiviral state. However, 
viruses have developed diverse strategies to inhibit type I IFN signaling. STAT2, an 
essential component of ISGF3 activity, is a common target hijacked by viral proteins. 
For example, the NS5 protein of flaviviruses, such as DENV and ZIKV, target 
STAT2 for degradation (30, 31). The HCV core protein was also reported to interact 
with STAT2 to disturb the function of ISGF3 complex (32). Phosphorylation of 
STAT2 was inhibited by C protein of Sendai virus, and yellow fever virus NS5S 
protein could prevent STAT2 binding to ISRE promoter (38, 39). However, the 
relationship between proteins encoded by CoVs and IFN signaling inhibition remains 
largely unclear. Limited research has shown that SARS-CoV nsp1 and ORF6 target 
STATI through inhibition of its phosphorylation or nuclear translocation, respectively 
(40, 41). SARS-CoV PLpro can disrupt the interaction between ERK1 and STATI to 
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antagonize type I IFN signaling (42). Here, we established that PDCoV nsp5, a 
3C-like protease, targeted STAT2 for cleavage. Furthermore, our study demonstrated 
that PDCoV nsp5 disrupted type I IFN signaling by cleaving STAT2, which was 
contributed by PDCoV nsp5 protease activity, independent of the ubiquitin 
proteasome pathway or apoptotic pathway. Importantly, endogenous STAT2 in 
PDCoV-infected cells was cleaved, which is consistent with the effect of nsp5 
overexpression. These results confirm the antagonistic property of PDCoV nsp5 in 
type I IFN signaling, which may affect the IFN response after PDCoV infection. 

We previously demonstrated that PDCoV nsp5 cleaves NEMO, a vital adaptor 
of the IKK complex, to reduce type I IFN production (33). The role of PDCoV nsp5 
in both IFN production and IFN signaling blockage suggests that PDCoV nsp5 limits 
the IFN-induced antiviral state in infected cells and contributes to rapid spread of 
PDCoV from an infected cell to neighboring cells. However, the relative impact of 
either NEMO or STAT2 in cells with PDCoV infection is unclear. Based on the 
identified residues of STAT2 and NEMO cleaved by PDCoV nsp5, to establish a cell 
line that expresses forms of NEMO and STAT2 with mutated cleavage sites will be 
helpful to evaluate the functional significance of potential cleavage of STAT2 or 
NEMO during PDCoV infection. In addition, NEMO and STAT2 are vital molecules 
for IFN production and IFN signaling, respectively. Previous studies have 
demonstrated that NEMO can be cleaved by the 3C or 3C-like proteases of PEDV, 
Hepatitis A virus (HAV), Food-and-mouth disease virus (FMDV), and Porcine 
reproductive and respiratory syndrome virus (PRRSV) (43-46). To investigate 
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whether STAT2 cleavage is a shared attribute among 3C-like proteases of different 
coronaviruses, we analyzed the nsp5 of two porcine coronaviruses, PEDV and 
transmissible gastroenteritis virus (TGEV), and five human coronaviruses (HCoVs) 
including SARS-CoV, Middle East respiratory syndrome virus (MERS-CoV), 
HCoV-229E, HCoV-OC43, and HCoV-NL63. We found that porcine STAT2 could not 
be cleaved by the nsp5 of PEDV and TGEV, and human STAT2 could not be cleaved 
by the nsp5 of all tested human coronaviruses (data not shown). However, whether 
cleavage of STAT2 is unique to PDCoV nsp5 requires further evidence in other 
coronaviruses. 

It is well known that STAT2 exhibits several functionally conserved regions, 
such as the critical tyrosine residue (Y690), which is phosphorylated to mediate STAT 
heterodimerization and transcription activation, and the transcriptional activation 
domain (TAD) at the C-terminal end, which is responsible for transcriptional activity 
after the ISGF3 complex binds to DNA (47-49). Our present study identified that 
PDCoV nsp5 cleaves STAT2 residues Q685 and Q758. It is considered that cleavage 
of Q685 by PDCoV nsp5 could disrupt the function of Y690 and inhibit STAT2 
activation. This finding also helps explain impaired tyrosine phosphorylation of 
STAT2 under PDCoV nsp5 overexpression upon IFN-a stimulation. The second 
cleavage of Q758, which is located in the TAD of STAT2, is inevitable to further 
disrupt the TAD, which helps to recruit other transcription co-factors. This may help 
explain why other cleavage products or the mixture of cleaved products of STAT2 
along with STAT! and IRF9 failed to efficiently activate ISRE promoter activity. 
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Taken together, cleavage at Q685 and Q758 by PDCoV nsp5 damaged the structural 
and functional integrity of STAT2, and affected subsequent STAT2 phosphorylation 
and ISGs induction. 

According to previous studies on CoVs, the general substrate recognized by 
CoV nsp5 was determined to be mainly the Pl, P2, and PI’ residues (50, 51). As 
shown by the analysis of PDCoV nsp5 substrate specificity (Fig. 5A), the P1 residue 
was a conserved Q residue and the P2 residue was primarily a leucine (L), which has 
a large side chain. Additionally, small residues generally occupy the P1’ position. The 
peptides EKANLQ¢gsJERKKYL and LEPMLQ7s3|STLEPV (P6 to P6’'; numbering 
based on STAT2) recognized by PDCoV nspS at the P1 and P2 positions are generally 
conserved substrates for all CoVs. However, glutamic acid (E) at the P1’ position, 
which shared a large negatively-charged side chain responsible for Q685 cleavage, 
was unusual. Therefore, the homology model for PDCoV nsp5 in complex with these 
two peptide substrates was constructed based on the structure of PEDV nsp5 (PDB 
identifier 4ZUH) to further investigate the changes at the P1’ position between the two 
substrates of STAT2. The PDCoV nsp5 S1’ subsite is composed of S25, A26, L27, 
C38, P39, H41, and V42 (Fig. 7A and 7B). The model showed that the long side chain 
of P1'-E fit comfortably in the pocket. Detailed molecular interactions between the 
two substrates and PDCoV nsp5 showed that the conserved Q at the P1 position was 
stabilized by two hydrogen bonds with H162 and E165, and the main chain of 
P2’-residue formed two hydrogen bonds with the main chain of A26 (Fig. 7C and 7D). 
Previous studies showed that PEDV nsp5S tolerates the positively-charged side chain 
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of arginine (R) of PEDV N protein and a larger hydrophobic residue valine (V) of 
NEMO at the P1’ position for cleavage (36, 43, 52). Our analysis revealed that the S1’ 
subsite of PDCoV nsp5 demonstrates flexibility to tolerate residues with a large 
negatively-charged side chain instead of small aliphatic residues, suggesting a broader 
selection in PDCoV nsp5 substrate recognition. Thus, the complete crystal structure of 
nsp5 and the cleaved peptide substrates of STAT2 could provide more direct evidence 
and detailed information of molecular interactions between the substrate and nspS, 
aiding in the design of targeted drugs against PDCoV nsp5. 

Similar to PDCoV, PEDV is a swine enteropathogenic CoV. Of 23 
PEDV-encoded proteins, at least 10 proteins have been identified as IFN antagonists 
(53). The mechanisms utilized by PEDV nspl, nsp5, PLP2, and N protein to 
antagonize type I IFN production have been elucidated (43, 53-55). As a newly 
emerging enteric CoV, the molecular mechanism by which PDCoV regulates IFN 
activity remains largely unknown. To date, only PDCoV nsp5 has been identified as 
an antagonist of both type I IFN production (33) and type I IFN signaling (present 
study). The further identification of PDCoV-encoded IFN antagonists and 
understanding of the interaction between PDCoV and the IFN signaling pathway will 


help to develop novel therapeutic targets and more effective vaccines. 


MATERIALS AND METHODS 
Cells and viruses. HEK-293T and PK-15 cells, obtained from the China Center for 
Type Culture Collection, were cultured at 37°C in 5% COz in Dulbecco's Modified 
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Eagle's medium (Invitrogen, USA) supplemented with 10% fetal bovine serum. 
LLC-PK1 cells for PDCoV infection were acquired from the ATCC (ATCC number 
CL-101) and cultured at the same condition described above. PDCoV strain 
CHN-HN-2014 (GenBank accession number KT336560) was isolated from a piglet 


with severe diarrhea in China in 2014 (56). 


Plasmids and luciferase reporter assay. Nsp5 gene of PDCoV strain CHN-HN-2014 
was amplified and cloned into pCAGGS-HA-C with a C-terminal HA tag. PDCoV 
nsp5 mutants were also constructed in the pCAGGS-HA-C background. Porcine 
STATI, STAT2, IRF9, JAK1, and TYK2 were cloned into pCAGGS-Flag with an 
N-terminal Flag tag. All STAT2 substitution mutants were constructed in the 
pCAGGS-Flag background by overlap extension PCR using specific mutagenic 
primers (available upon request). STAT2 deletion mutants were also cloned into 
pCAGGS-Flag. All constructed plasmids were confirmed by sequencing. 

The luciferase reporter plasmid ISRE-Luc used in this study has been previously 
described (57). For luciferase reporter assays, HEK-293T cells were transfected with 
reporter plasmid ISRE-Luc and pRL-TK. Twenty-four hours after transfection, the 
cells were stimulated with IFN-a (Catalog no. 11101-2, PBL Assay Science) for 12 h 
at a final concentration of 1,000 U/ml. The lysed cells were prepared and the 
luciferase reporter assay system (Promega, Madison, WI) was utilized to determine 
the luciferase activities of lysed cells. The activities were normalized to that of 
pRL-TK (Promega). 
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RNA extraction and quantitative real-time PCR. To determine the effect of 
PDCoV nsp5S on the expression of ISGs, including ISG15, ISG54, ISG56, and OAS1, 
HEK-293T cells in 24-well plates were transfected with 1 ug of empty vector or nsp5 
expression plasmid. After 24 h, the cells were treated with IFN-o (1,000 U/ml) for 12 
h. To evaluate the ability of the cleavage products of STAT2 to induce ISGs, 
full-length STAT2 or its deletion mutants were cotransfected with STAT1 and IRF9 
into PK-15 cells for 32 h. Total RNA was extracted from transfected cells with TRIzol 
reagent (Invitrogen). RNA was then reverse transcribed into cDNA by avian 
myeloblastosis virus reverse transcriptase (TaKaRa, Japan). Quantitative real-time 
PCR (qPCR) experiments were performed in triplicate. mRNA expression levels were 
normalized to that of glyceraldehyde-3-phosphate dehydrogenase (GAPDH). All 


qPCR primers used in this study are listed in Table 1. 


Western blot analyses. Cells were cultured in 60-mm dishes and harvested with lysis 
buffer (Beyotime, China). The samples were then resolved by SDS-PAGE and 
transferred to polyvinylidene difluoride membranes (Millipore, USA) to determine 
protein expression. Cells were also treated with a broad caspase inhibitor 
(Z-VAD-FMK; Beyotime, China) or a proteasome inhibitor (MG132; Beyotime, 
China) at a final concentration of 20 uM. The overexpression of JAK1, TYK2, IRF9, 
STATI, STAT2, and STAT2 mutants was evaluated using an anti-Flag antibody 
(Macgene, China). An anti-HA antibody (MBL, Japan) was used to analyze the 
expression of PDCoV nsp5. STATI polyclonal antibody (Santa Cruz, USA), STAT2 
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antibody (ABclonal, China), JAK1 antibody (Cell Signaling Technology, USA), and 
TYK2 antibody (Cell Signaling Technology, USA) were utilized to detect each 
respective endogenous protein. Phospho-STAT1 antibody (Cell Signaling Technology, 
USA), Phospho-STAT2 antibody (ABclonal, China), Phospho-JAK1 antibody (Cell 
Signaling Technology, USA), and Phospho-TYK2 antibody (Cell Signaling 
Technology, USA) were utilized to identify the phosphorylated form of each 


respective endogenous protein. 


Homology modeling. A BLASTP search of PDCoV nsp5 sequence against the 
Protein Data Bank (PDB) library revealed PEDV nsp5 (PDB identifier 4ZUH) to be a 
close homolog. Thus, the homology model of PDCoV nsp5 was generated based on 
the crystal structures of PEDV nsp5 using MODELLER 9.14 software (58). The 
distribution of hydrophobic and hydrophilic interactions between the substrates and 


protease was created by the LigPlot+ program. 


Statistical analysis. All experiments were performed in triplicate. Significant 
differences were determined using Student’s ¢ test. P values <0.05 were considered 


Statistically significant. 
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FIGURE LEGENDS 

FIG 1. PDCoV nsp5 antagonizes type I IFN signaling. (A) HEK-293T cells 
cultured in 24-well plates were transfected with PDCoV nsp5 expression plasmid or 
vector. After 24 h, cells were treated with 1,000 U/ml of IFN-a for 8 h and analyzed 
by qPCR. (B, C and D) HEK-293T cells (B), PK-15 cells (C) or ST cells (D) cultured 
in 24-well plates were transfected with various concentrations of PDCoV nsp5 
expression plasmid along with ISRE-Luc plasmid and pRL-TK plasmid. After 24 h, 


cells were treated with 1,000 U/ml of IFN-a for 12 h, followed by luciferase assays. 


FIG 2. PDCoV nsp5-mediated inhibition of type I IFN signaling requires its 
protease activity. (A) Sequence alignment of nsp5 for the four genera of CoVs. 
Delta-CoVs were represented by PDCoV (accession number ALS54085.1) and 
HUK13 (accession number YP_002308505.1); alpha-CoVs were represented by 
PEDV (accession number NC_004004) and NL63 (accession number AFV53147.1); 
beta-CoVs were represented by SARS-CoV (accession number NP_828850.1) and 
MHV (accession number NP_068668.2); and gamma-CoVs were represented by IBV 
(accession number NP_740623.1) and TCoV (accession number YP_001941178.1). 
Secondary structures of CoV nsp5s were analyzed using the ESPript website 
(http://espript.ibcp.fr/ESPript/ESPript/index.php). (B) Structural alignment of the 
conserved residues His41 and Cys144 in CoV nsp5s. The homology model of PDCoV 
nsp5 was generated based on the crystal structures of PEDV nsp5 by MODELLER 
9.14 software. The 3D structures of PEDV nsp5 (PDB identifier 4ZUH), SARS-CoV 
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nsp5 (PDB identifier 2Q6G), and IBV nsp5 (PDB identifier 2Q6D) were obtained 
from the Protein Data Bank. Red arrows indicate conserved enzymatic proteolysis 
residues. (C) HEK-293T cells in 24-well plates were transfected with PDCoV nsp5 
expression plasmid or its protease-defective mutants C144A and H41A to analyze 


ISRE promoter activity, as described in Fig. 1B. 


FIG 3. PDCoV nsp5 inhibits type I IFN signaling by targeting STAT2. (A) 
PDCoV nsp5 does not affect JAK1 and TYK2 phosphorylation and expression. PK-15 
cells were transfected with PDCoV nsp5 expression plasmid or vector. After 24 h, 
cells were treated with IFN-a (1,000 U/ml) for 1 h and collected for western blot 
analysis. JAK1 antibody, TYK2 antibody, Phospho-JAK1 (P-JAK1) antibody and 
Phospho-TYK2 (P-TYK2) antibody were utilized to detect each respective 
endogenous protein. (B) The level of endogenous STAT2 and phosphorylation were 
both reduced under the expression of PDCoV nsp5. PK-15 cells were transfected with 
PDCoV nsp5 expression plasmid or vector and treated as described in panel A. STAT1 
polyclonal antibody, STAT2 antibody, Phospho-STAT1 (P-STAT1) antibody, 
Phospho-STAT2 (P-STAT1) antibody were used for western blot. (C) HEK-293T cells 
were cotransfected with Flag-appended JAK1, TYK2, STAT1, STAT2, or IRF9 
expression plasmids along with PDCoV nsp5. Cells were lysed at 32 h after 


transfection and analyzed by western blotting. 


FIG 4. PDCoV nsp5 cleaves STAT2 to disrupt type I IFN signaling. (A) 
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HEK-293T cells were cotransfected with various concentrations of PDCoV nsp5 
expression plasmid and Flag-appended STAT2. Cells were then lysed for western blot 
analysis at 32 h after transfection. (B) HEK-293T cells were cotransfected with 
wild-type PDCoV nsp5 or its protease-defective mutants (C144A and H41A) and 
STAT2 expression plasmid as described in panel A. (C) HEK-293T cells were 
cotransfected with wild-type PDCoV nsp5 and STAT2 expression plasmids. After 16 
h, cells were treated with MG132 or Z-VAD-FMK (final concentration of 20 4M) for 
8 h. Cell lysates were prepared and analyzed by western blotting. (D) LLC-PK1 cells 
were infected with PDCoV strain CHN-HN-2014 (MOI=2) and then harvested for 
STAT2 mRNA detection by qPCR at 18 h post-infection. (E) LLC-PK1 cells were 
infected with PDCoV as described in panel D and then lysed for western blot analysis 


at 18 h post-infection. 


FIG 5. Q685 and Q758 are the sites of PDCoV nsp5-mediated STAT2 cleavage. 
(A) Sequence logo of the polyprotein junctions cleaved by PDCoV nsp5. An amino 
acid sequence logo of the substrate was generated by WebLogo 3 
(http://weblogo.threeplusone.com/). (B) The estimated STAT2 cleavage product by 
PDCoV nsp5. HEK-293T cells were transfected with STAT2 (1-637 aa), STAT2 
(1-719 aa), or STAT2 (1-758 aa) deletion mutant. Cells were then harvested and 
analyzed by western blot analysis. (C) Schematic representation of STAT2 and its 
mutants. (D) HEK-293T cells were transfected with PDCoV nsp5 and wild-type 
STAT2 or STAT2 mutants as indicated. Cells were lysed at 32 h after transfection and 
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analyzed by western blotting. (E) HEK-293T cells were transfected with PDCoV nsp5 


and wild-type STAT2 or mutant STAT2-Q685A-Q758A as described in panel D. 


FIG 6. PDCoV nsp5-mediated STAT2 cleavage is involved in inhibition of type I 
IFN signaling. (A) PK-15 cells were cotransfected with STAT1, IRF9, and wild-type 
STAT2 or the indicated deletion mutants, along with ISRE-Luc plasmid and pRL-TK 
plasmid. After 32 h, cells were harvested for luciferase assays. (B, C and D) PK-15 
cells were transfected with expression plasmids of STAT1, STAT2, and IRF9 as 
described in panel A. After 32 h, cells were harvested and the mRNA levels of JSG/5 


(B), ISG54 (C) and OAS/ (D) were detected by qPCR. 


FIG 7. Homologous modeling of PDCoV nsp5 with the two cleaved peptide 
substrates. (A and B) Homology model of PDCoV nsp5 in complex with the cleaved 
peptide substrates (STAT2) based on the structure of PEDV nsp5 (PDB identifier 
4ZUH). The residues EKANLQ¢ssERKKYL (A) or LEPMLQ7;sSTLEPV (B) (P6 to 
P6’; numbering based on STAT2) occupy, and thereby define, the subsites S6 to S6’ of 
PDCoV nsp5. The positions of P6 to P6’ and potential residues forming S1’ (red) 
subsites are labeled. (C and D) The distribution of hydrophobic and hydrophilic 
interactions between PDCoV nsp5 and EKANLQ¢gsERKKYL (C) or 
LEPMLQ7ssSTLEPV (D). Carbon, nitrogen, and oxygen atoms are marked as black, 
blue, and red circles, respectively. Hydrogen bonds are shown as green dashed lines 
labeled with the distance between the donor and corresponding acceptor atom. Arcs 
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Table 1. Primers used for qPCR 


Primer Sequences (5' to 3') 
h-ISG15-F GGGACCTGACGGTGAAGATG 
h-ISG15-R CGCCGATCTTCTGGGTGAT 
h-ISGS54-F CACCTCTGGACTGGCAATAGC 
h-ISG54-R GTCAGGATTCAGCCGAATGG 
h-ISGS6-F GCTTTCAAATCCCTTCCGCTAT 
h-ISG5S6-R GCCTTGGCCCGTTCATAAT 
h-OAS1-F CGTGTTTCCGCATGCAAATC 
h-OAS1-F GCGAACTCAGTACGAAGCTG 
h-GAPDH-F TCATGACCACAGTCCATGCC 
h-GAPDH-R GGATGACCTTGCCCACAGCC 
h-B-actin-F CCTGGCACCCAGCACAAT 
h-B-actin-R GCTGATCCACATCTGCTGGAA 
P-STAT2-F TGGATGGAGAGTTGGAGCAG 
P-STAT2-R CTTGCTCCCAGTCTTGAGGA 
p-ISGI15-F CCTGTTGATGGTGCAAAGCT 
p-ISG15-R TGCACATAGGCTTGAGGTCA 
p-ISGS54-F CTGGCAAAGAGCCCTAAGGA 
p-ISG54-R CTCAGAGGGTCAATGGAATTCC 
p-OAS1-F AAGCATCAGAAGCTTTGCATCTT 
p-OAS1-R CAGGCCTGGGTTTCTTGAGTT 
p-GAPDH-F ACATGGCCTCCAAGGAGTAAGA 
p-GAPDH-R GATCGAGTTGGGGCTGTGACT 
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Substrate: Leu1(P6)-Glu2(P5)-Pro3(P4)-Met4(P3)-LeuS(P2)- 
Gin6(P1)(758)-Ser7(P1’)-Thr8(P2’)-Leu9(P3’)- 
Glu10(P4’)-Pro11(P5’)-Val12(P6’) 
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